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Abstract 

We study the influence of fluctuations in molecular shape on the stability of the biaxial nematic 
phase by generalizing the mean field model of Mulder and Ruijgrok [Physica A 113, 145 (1982)]. 
We limit ourselves to the case when the molecular shape anisotropy, represented by the alignment 
tensor, is a random variable of an annealed type. A prototype of such behavior can be found 
in lyotropic systems - a mixture of potassium laurate, 1-decanol, and D2O, where distribution of 
the micellar shape adjusts to actual equilibrium conditions. Further examples of materials with 
the biaxial nematic phase, where molecular shape is subject to fluctuations, are thermotropic 
materials composed of flexible trimeric- or tetrapod-like molecular units. Our calculations show 
that the Gaussian equilibrium distribution of the variables describing molecular shape (dispersion 
force) anisotropy gives rise to new classes of the phase diagrams, absent in the original model. 
Depending on properties of the shape fluctuations, the stability of the biaxial nematic phase can 
be either enhanced or depressed, relative to the uniaxial nematic phases. In the former case the 
splitting of the Landau point into two triple points with a direct phase transition line from isotropic 
to biaxial phase is observed. 

PACS numbers: 61.30.Cz, 64.70.Md, 05.70.-a 



*e-mail addressponga@th.i f.uj.edu.pl 
'I' e-mail addressjwydro@univ-metz.fr 



1 



I. INTRODUCTION 



The biaxial nematic phase is one of the perennially challenging problems of experimental 
soft matter physics. Although predicted theoretically by Freiser in 1970 l| it was not 
until 1980 that the first experimental observation of this phase has been reported The 
phase was discovered in a lyotropic, ternary mixture of potassium laurate, 1-decanol, and 
D 2 and its stability attributed to shape change of the micellar aggregates as function of 
temperature and concentration of amphiphilic molecules 3J. The search for a thermotropic 
biaxial nematic has proved highly controversial for more than two decades. A first qualitative 
report on a stable uniaxial-to-biaxial nematic phase transition has been reported by Li et. 
al. {4] in their system of flexible, ring-shaped trimeric liquid crystal. Recently, this phase 
has also been detected in two classes of thermotropic materials. It was reported in 'banana- 
shaped' mesogens 

BB 

and in liquid crystalline tetrapodes [7j. The last class of materials 
is particularly interesting for it comprises molecules with four mesogenic units connected 
together through a flexible spacer. The optimal packing of such tetrapodes in the biaxial 
nematic phase is achieved for a quasiflat, platelet configuration of the tetrapod's mesogens 
that are, on the average, tilted in the plane of the platelet. 



A challenge for theory is to find molecular factors responsible for absolute stability o 
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the 
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observed biaxial nematic phase. The presently existing microscopic models [l, 9 
show that the molecular shape- and pair interaction biaxiality are obviously important for 
;he biaxial phase to exist. However, as numerous experimental reports have demonstrated 



12j, they seem not sufficient to get the absolutely stable biaxial phase. In the present paper 



we show that a variable (fluctuating) anisotropy of the molecular shape can be an important 
stabilizing factor for this phase to occur. Indeed, it seems that this factor appears commonly 
in the micellar- trimeric- and tetrapod systems. Let us mention that the theoretical studies 
and computer simulations so far have addressed the size and shape of the micelles in lyotropic 

systems 



but a connection between the shape anisotropy distribution and the 



stability of the biaxial nematic phase have not yet been explored. 

The present paper is arranged as follows. After introduction of the model in Sec. II, we 
explore stability of the biaxial nematic phase on shape fluctuations in Sec. III. Section IV 
is devoted to a short summary. 
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II. THE MODEL 



We assume that the hamiltonian if ({/i}, {S}) of N liquid crystalline molecules with 
dynamical degrees of freedom {S} also depends on the {fi} variables, which parameterize 
molecular shape. The partition function is calculated for each allowed configuration of 
randomly chosen {/i}. Then, in analogy to statistics of disordered systems with annealed 
disorder [l7|. {/i} is treated as a set of dynamical variables that adjust to actual equilibrium 
conditions. Under these circumstances the free energy, F, of the system can be approximated 
by the logarithm of the {/i} averaged partition function 



where P({fi}) is the probability distribution of the {//} variables. A role played by the 
distribution P on the formation of the biaxial nematic phase is studied by generalizing a 
very elegant mean- field model of Mulder and Ruijgrok (MR), which employs a connection 
between the SU(3) symmetry and the biaxial nematic liquid. The most important feature 
of the model is that its partition function can be calculated exactly, which, as we are go- 
ing to show, allows one for a semi-analytical treatment of the annealed average ([2]). More 
specifically, in the MR model the internal, dynamical state of each molecule is parameter- 
ized by eight degrees of freedom: three components L a of the angular momentum L and 
five components Q a p of the symmetric and traceless quadrupole moment Q. These eight 
variables are next combined to form eight generators S a of the SU(3) algebra, establishing 
equivalence between {L, Q} and S = Ylt=i S a ^ a , where A a are the Gell-Mann matrices: 
(S) a /3 = Q a f3 — f J2-y £ ai3- y L 1 , a, (3, 7 = 1 . . . 3. The MR Hamiltonian is the mean-field (MF) 
version of the Heisenberg-type interaction 



F 




—knT In Z 



(1) 
(2) 




(3) 



i=l 
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where the dynamical variables Sj are subject to two i-independent constrains, represented 
by Casimir invariants of the SU(3) algebra 

Tr(S 2 ) = Tr(^) = ][>J Q = 2/ 2 ( Mi ) (4) 



a=l 
3 



Tr(Sf ) = Tr(Q, 3 ) = £ ^ a = 2h(jn) (5) 

a=l 

with fj,i i + yU ij2 + Hi,3 — 0; N is the number of molecules and S is the thermodynamic average 
of Sj. The /i^Q, variables are the eigenvalues of the traceless matrices Sj, or Qj, where Qj 
is the the quadrupole moment of the i — th molecule at rest (Lj=0), obtained from Sj by 
applying an SU(3) transformation. Hence, if Qj represents e.g. the quadrupole moment of a 
mass distribution the intrinsic properties of a given molecule, like the ratios of its principal 
axes, enter through the Casimir invariants ^(fii), ^(/Ltj), Eqs. (T4151) . This remains in full 
analogy with what we practise for an ordinary quadrupolar tensor, where the Lr and / 3 -like 
invariants are used to ch— biaXlallty of levant p hysi c al observe QQ. Mo r e 
specifically, depending on the values of I2 and ^3, or their ratio 



v^-^t (6) 



w 

^ )2 



three possibilities can be distinguished: (a) for I2 = I3 = the tensor is isotropic; (b) for 
3/f = /f the tensor is uniaxial and (c) for 3/f < /f the tensor is biaxial with maximal 
biaxiality being obtained for ^3 = 0. The sign of ^3 decides about whether the tensor is 
prolate (plus sign) or oblate (minus sign). Respectively, w = 1 (w = — 1) and \w\ < 1 
refers to rod-like (disk-like) and biaxial molecules. By construction the MR model is SU{3) 
invariant with degrees of freedom running over the group manifold and its free energy is 
given in an analytical form as derived by Itzykson and Zuber j^ol ]. 

From physical point of view the model matches the standard Landau-deGennes phase di- 



agram for biaxial nematics [18j, known to characterize systems with purely dispersion-type 
of interactions. Jt^also reproduces the mean- field results for the dispersion model of Luck- 
hurst et al. 
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22], Fig. [TJ given we take the pair interactions of the form V = — eTr(RR'), 
where R denotes the normalized quadrupole tensor (TrR 2 = 1) defined through the rela- 
tion: \/T+ 2/t 2 R = ^g(3z <g> z — 1) ± k(x <S> x — y <g> y); k is the ratio of biaxial-to-uniaxial 
polarizability of the molecule. The relative error for the Njj — Nb boundary calculated for 
both models does not exceed 2% and is even smaller for the nematic order parameter (see 



Fig. 4 in |8|). Note however the usefulness of the SU(3) symmetry. It avoids calculation 
of integrals over Euler angles inherent to the dispersion model, which, in turn, allows for 
comprehensive studies of flexible quadrupoles. 

A generalization of the MR model to systems with variable anisotropy of the molecular 
shape (dispersion forces) is straightforward. We assume that /i^x and /ij j2 are annealed de- 
grees of freedom weighted with the probability PQu^i, fj,i >2 , 112,1, /Mr, 2) ~ YliLi P{lH,i, Hip)- 
The mean-field partition function, Z = Z± , and the dimensionless free energy are then given 
by 

Tr(SS) 

Zi = e~— Q (7) 
Q = (JdS e^Si^fl - h(^))S(^p - J 8 (/i))> (8) 
F = -Nt logZi, (9) 

where t = k B T / J is the dimensionless temperature and where (...) = J (...)P(/i 1; /x 2 )d/xi<i/X2- 
According to Itzykson and Zuber [2^] the integral over S in (jSJ) can be carried out to give 

where 

/ 



A 7 W 



D = det 



e 7lMl g72Ml e 73Ml 
e 7lM2 g72M2 e 73M2 ; (11) 

y e^ 1M3 e^ 2A13 e^ 3M3 y 



and where 

A x = (xi -x 2 )(x 2 -x 3 )(x 3 -zi). (12) 

The eigenvalues % of the traceless matrix S/(2t), are determined from the minimum of 
the free energy (jHJ). In analogy to ( T4l5i) the invariants /2C7) an d ^3(7) are used to distin- 
guish between (a) the isotropic phase (^2(7) = ^3(7) — 0); (b) the uniaxial nematic phase 
(3/ 3 (7) 2 = ^(7) 3 ) an d ( c ) the biaxial nematic phase (3/ 3 (7) 2 < /2(7) 3 )- In addition, 13(7) 
is positive for prolate uniaxial phase (Nu+) and negative for oblate uniaxial phase (iV7j_). 

The annealed averaging (...) over P(//i,/i 2 ), Eq. (|S]), has a very simple interpretation 
in the mean-field theory. Setting %, = 0, (a = 1, 2), which is always one of the mean-field 
solutions, we find that P(/ti, ^2) matches the density distribution of the molecular shape 
anisotropy in the reference (stable or metastable) disordered phase. We believe therefore 
that for a credible choice of P(/ii, ^2) the model correctly reproduces generic phase behavior 



for flexible quadrupoles in the vicinity of the isotropic phase. Clearly, the original MR model 
is recovered if P(/xi,/X2) is given by Dirac delta distribution. In what follows we take P to 
be the Gaussian distribution. This choice is consistent with the maximum entropy principle 
applied in the isotropic phase and the observation that usually only first two moments of 
P can be determined reasonably well from experiment {j^. Assuming that in the reference 
(disordered) phase these moments are given by: < \i a >= m a , < (fi a — m a ) 2 >= a 2 a and 
< (fi\ — rrii)(fi2 — vn-i) >= Action (a = 1, 2) we find 

P(fj, u /i 2 ) = V°L e-^^^-^K^-^), (13) 

where 



(71<T2 



with A being the real parameter, a=|(l + \J\ + 4A 2 ) and — 1 < A=^ < 1- The distribution 
P(fJ>3) of /i3, obeying the constrain + fi 2 + A*3 = 0, is also a Gaussian with average 
< fi^ >= m 3 = —mi — and dispersion cr| = o\ + 2Acri(T2 + of — cr 2 1 < 0"3 < 
°"i + °"2) : P{^z) — ^Kf- ex P(~ ^±Z^+m) ), i n general, the parameters of the distribution 
( TT3"j) can depend on temperature, but this possibility, which can be relevant for quantitative 
understanding of phase diagrams in lyotropic systems, will not be discussed here. 

III. RESULTS 

The phase diagrams, obtained from the global minimization of the free energy with 
respect to {7^, a = 1,2}, depend on the values of five parameters: mi, 7712, c 2 and A (or 
03) of the Gaussian distribution ([TBI . To make a direct comparison with the earlier work 8j 
we use, instead of mi and m.2, the molecular shape parameter w(m) and ^(m), Eqs. (I4f6j) . 
A connection between the parameterizations is given by 

mi + m 2 = u (15) 
mim 2 = u 2 - J 2 , (16) 

where u is a solution of the cubic equation 

u 3 -u/ 2 + ^(/ 2 ) 3 /\; = 0. (17) 
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The three real roots U\ < w 2 < W3 of Eq. (fTTj) correspond to different permutations between 
axes of the molecule-fixed frame. In what follows we choose the ^-solution to identify the 
corresponding mi and m 2 . Solving Eqs. (11511161) for given w 2 and 7 2 still leaves a freedom to 
choose the pair {m l5 m 2 } up to a permutation. We select the solution for which mi < m 2 . 
Permutation symmetry of Q: Q({m a }, {a a }, ...) = Q({mp^ a )}, {ap^ a )}, ...), where V is an 
arbitrary permutation of {1, 2, 3}, allows us to construct the remaining diagrams for m\ > m 2 
from the ones given. 

Numerical calculations are carried out for fixed values of <Ti, <t 2 and A. The phase diagrams 
are shown in plane of the molecular shape parameter w(m) and the reduced temperature 
t/ti, with ti > </2 g M ^ > being the isotropic-nematic transition temperature for w{m) = 0. 
In all cases the numerical value of J 2 was fixed to 7 2 = \- The diagram for J 2 = |£ 2 , <Tj and 
A, where £ is an arbitrary real number, can be obtained from that for I 2 — j, Oi — > crj/^ and 
A, which follows from invariance of Q, Eq. ([8]), with respect to £-rescaling of the parameters: 
Q({rrii}, {di}, \,t) = Q({^,mi}, {£(7j}, A, £t). In addition, the invariance of Q, Eq.Q, with 
respect to change of {/j,, 7} into {— /x, —7} makes the phase diagrams symmetric with respect 
to the line w(m) = 0. 

Numerical minimization of the free energy allows us to divide all the diagrams into classes 
shown in Figs. [20 The corresponding isotropic-nematic transition temperature tL is plotted 
in Figs. I5ll6l At high temperatures, in the vicinity of isotropic-nematic phase transition, the 
results can be understood qualitatively from the expansion of Q, Eq. ( JTOl) . about 7^ = 0. Up 
to sixth-order in 7^ it reads 



Q « l + -(/ 2 )/ 2 ( 7 ) + T7 :(/3)/3(7) + ^(/2)/2(7) 2 



+ —(hh) h(l) h(l) + f (Io) — (A 2 ) 1 hil? 

?0 \ a 1 2\jj V 12 096 V 2 1 5040 V 3 7 V /; 




(18) 
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where the averages over fi are given by 



(h) — I20 + ^02 



( I 3 ) = / 30 + 3/ 12 
{ID = J 2 2 + 6J 22 + 3/, 



04 



( h h ) = hoho + 6J20I12 + 4/30/02 + 9/u 



(19) 



( II ) = 4) + 184 + y^oe + 15/i 2 (/ 30 + 3/ 12 ) 




with 



/, 



pg 



(/ 2 = / 20 , / 3 = /30). 



(20) 



Using Eqs. (|9]fT8T) we find that in the limit of small {%} the bifurcation from isotropic 
to nematic phase takes place at tf,=|(/ 2 ). Due to the /o 2 contribution to ( / 2 ) (/o 2 > 
0) the bifurcation temperature, tb, is always greater than the corresponding bifurcation 
temperature for the mono-dispersive system with molecules characterized by the average 
shape parameters {m,}. Similar result should hold for the transition temperatures, for they 
usually follow behavior of 

In the vicinity of the isotropic phase the terms higher than sixth order in {7^} can be 
neglected in the expansion (1181) . Out of the six terms that are left one can associate ( / 3 ) with 
an effective molecular shape anisotropy of the system in the isotropic phase. The lyi term, 
contributing to and being of undetermined sign, effectively changes this anisotropy and 
thus has a profound effect on stability of the biaxial phase. More specifically, as ( / 2 ) > / 2 o 
and (/|) > /| , for / 30 and /i 2 of opposite sign with |/ 12 | < |/ 30 | the ( / 3 ) coefficient can 
effectively be reduced by the shape fluctuations. If, simultaneously, the invariant ( / 2 /3 ) 
is also reduced by fluctuations, which as we checked is easily achieved in the parameter 
space, the range of stability of the biaxial phase in (w(m), t/ti) space gets enhanced as 
compared to the case without shape fluctuations. That is, a sufficient condition to observe 
a constructive role of shape fluctuations in stabilizing the biaxial nematic phase is the 
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simultaneous fulfilment of two inequalities 



\{hh)\ < ho\ho 



\(h)\ < |ii 



30 



(21) 
(22) 



Interestingly, any shape fluctuations about spherically symmetric shape (m a = 0) stabilize 
biaxial nematic phase of maximal biaxiality ((1^/3) = 0), without intermediate uniaxial 
phase. 

Now we turn to detailed analysis of the model. We carried out numerical minimization to 
determine phase diagrams for a representative set of model parameters. All distinct classes 
of the diagrams identified are shown in Figs. [20 In particular, we found that the class 
of parameters where the biaxial nematic phase enhances its stability is much more reacher 
that the condition ( 12T]) may suggest. However, we are unable to find numerical limitations 
on the model parameters in a compact form, except for some limiting cases. But even in 
these limiting cases we recover all observed classes of the phase diagrams. The simplest 
case occurs when one of the dispersions, say <7 2) vanishes. This corresponds to the Gaussian 
distribution for fi\ and fi^ (<7i = 03), and Dirac delta distribution for fi2- Phase diagrams, 
influenced by this polydispersivity, already exhibit quite different topology as compared to 
the original MR model. This we illustrated in Fig. [21 where diagrams B and C correspond 
to o"! = 0.15 and a± = 0.3, respectively. Note a considerable enhancement of stability of the 
biaxial phase along with splitting of the original MR quadruple Landau point into two triple 
points. The triple points are connected by the line of the first order transition between the 
isotropic and biaxial phases. Moreover, as expected, the transition temperature, t^, of the 
isotropic- nematic phase transition increases with increasing value of a±, Fig. [5j 

Subsequent case to consider is the full Gaussian distribution. We limit ourselves to the 
symmetric distributions with <j\ — a^. Now the phase diagrams exhibit yet another topology, 
which is represented by the diagram D in Fig. [31 Amazingly, there are two triple points and 
the Landau point on one diagram, and two lines of the direct first order phase transition 
between the isotropic and biaxial phases for |iu(m)| <J 0.8. Stable uniaxial phases form 
two bubble-like islands. As before, the transition temperature, ti, between the isotropic- 
and nematic phases is higher than that of the MR model and increases with increasing 
A, Fig. [6j Another possible diagrams obtained for this case are shown in Fig. [H In the 
diagram F the biaxial phase is practically eliminated being reduced to a line w(m) = 0. 
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It is interesting to follow reduction of stability of the biaxial phase by comparing diagrams 
B, H, G and F, where changes of 02 from zero to a\, for fixed <Ti, correspond to successive 
phase diagrams. Numerical estimates of low-temperature part of the phase diagrams have 
been checked to stay consistent with asymptotic expansion (| — > 00) of the selfconsistent 
equations {dF/d^i = 0} for {7i}. 

IV. SUMMARY 

We have studied the influence of the variable molecular shape anisotropy on stability of 
the biaxial nematic phase. To make the analysis as simple as possible we generalized the 
exact mean field solution obtained by Mulder and Ruijgrok for biaxial molecules to the 
case when the quadrupole tensor is a dynamical variable. We assumed that at equilibrium, 
the molecular shape anisotropy can be approximated by the (annealed) distribution, P, of 
the molecular parameters In the mean-field approximation the natural choice for 

P, consistent with the maximum entropy principle applied in the isotropic phase, is the 
two-dimensional Gaussian distribution. 

The nonzero second moments of the Gaussian distribution lead to important remodeling 
of the original MR phase diagram. We observe that polydispersivity changes the range of 
stable biaxial phase and that behavior of the system can qualitatively differ from its mono- 
dispersive counterpart characterized by the average shape parameters {mj}. Generally, 
the transition between the isotropic and the nematic phases occurs at higher temperatures 
when molecular shape changes are allowed. The phase diagram is modified, for instance, by 
showing the quadruple Landau point being splitted into two triple points connected by a line 
of first order transitions between the isotropic and biaxial phases. By comparing diagrams 
B, D and E we can conclude that strong correlations between shape fluctuations along 
main molecular axes (|A| — > 1), or fluctuations along two of the three molecular principal 
axes lead to particularly large region of stable biaxial phase. Importantly, fluctuations about 
isotropic shape give rise to stable biaxial nematic without intermediate uniaxial phase, while 
fluctuations about fixed w(m), diagram G, show on the temperature axis two uniaxial phases 
separated by the biaxial nematic. 

In some cases the biaxial phase can be destabilized in the vicinity of the isotropic phase 
giving only the uniaxial nematic phases and first-order phase transitions between them (class 
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F, Fig. |4]). Similar case has recently been observed by Bates 23| in his computer simula- 
tion of a generic, flexible V-shaped molecules on a lattice. The only difference between 



rat we do not observe a biaxial nematic phase at low 



23] . A reason for that is our neglecting of temperature 



our predictions and that of 23j is t 
temperatures, shown in Fig. 5(a) of 
dependence in the isotropic distribution (1T3"|) at low temperatures 24|. Clearly, to be con- 
sistent with general thermodynamics at T = the distribution (fl3|) should approach Dirac 
delta function about m a . In our studies we have disregarded any temperature dependence in 
(fT3|) . being primarily interested in system's behavior close to the isotropic phase. However, 
the diagram predicted in 23] can also be obtained within our model if we replace a a by 
\/ia a in (fTBl . Then the ground state of (J9j) for w(m) ^ ±1 would always be the biaxial 
nematic phase and, consequently, the phase diagrams of the class F, Fig. IH would go into 
generic diagrams found by Bates. 
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FIG. 1: Comparison of the MR phase diagram 
Luckhurst et al. dispersion model 



21 



(continuous lines) with one calculated for the 



22J (dashed lines); -l<w = \/6Tr(R 3 ) 



1 - fe2 . i < 1. 

(1+2k 2 )2 
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FIG. 2: Phase diagrams for 02 = 0. Diagrams labeled Ai, Bi and Cj correspond to o\ = 
0.15 and 0.3, respectively. Subscript i = 1 refers to the isotropic (Iso)-nematics phase transition 
lines whereas subscript i = 2 refers to the uniaxial nematic-biaxial nematic (Nb) lines. 
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FIG. 3: Phase diagrams for o\ = 02. Diagram labeled Ai corresponds to a\ = [8|], while diagrams 
Di and Ei to o\ = 0.15. Subscript i = 1 refers to the isotropic (Iso)-nematics phase transition lines 
whereas subscript i = 2 refers to the uniaxial nematic-biaxial nematic (Nb) lines. For diagrams 
Di and Ei we took A = 0.99 and A = —0.99, respectively. 
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FIG. 4: Phase diagrams for u\ = 0.15 and A = 0. Diagram B{ corresponds to 02 = 0, while 
diagrams Fi, Gi and Hi correspond to oi = 0.15, 0.1 and 0.05, respectively. Subscript i = 1 refers 
to the isotropic (Iso)-nematics phase transition lines whereas subscript i = 2 refers to the uniaxial 
nematic (Njj+, Njj— )-biaxial nematic lines. The lines of isotropic-nematic phase transitions occur 
in the following order (from left to right): B%, Gi, F\ and H\. 
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FIG. 5: Isotropic-nematic transition temperature tj_, for w(m) = as function of a\ or a 2, collec- 
tively denoted a. For continuous line a = a\ with 02 = whereas for dashed line a = 02, o\ = 0.15 
and A = 0. Points A, B and C correspond to a\ = 0, 0.15 and 0.3, respectively. Points F, G and 
H correspond to 02 = 0.15, 0.1 and 0.05, respectively. 
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FIG. 6: Isotropic-nematic transition temperature tj_, for w(m) = as function of A for o\ = 02- 
Points D, E and F correspond to A = 0.99, -0.99 and 0, respectively. 
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